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Abstract 



We are interested in a stochastic model of trait and age-structured population undergo- 
ing mutation and selection. We start with a continuous time, discrete individual-centered 
. population process. Taking the large population and rare mutations limits under a well- 

chosen time-scale separation condition, we obtain a jump process that generalizes the Trait 
Substitution Sequence process describing Adaptive Dynamics for populations without age 
' structure. Under the additional assumption of small mutations, we derive an age-dependent 

p I , ordinary differential equation that extends the Canonical Equation. These evolutionary ap- 

proximations have never been introduced to our knowledge. They are based on ecological 
phenomena represented by PDEs that generalize the Gurtin-McCamy equation in Demogra- 
phy. Another particularity is that they involve a fitness function, describing the probability of 
invasion of the resident population by the mutant one, that can not always be computed ex- 
plicitly. Examples illustrate how adding an age-structure enrich the modelling of structured 
population by including life history features such as senescence. In the cases considered, we 
^ ' establish the evolutionary approximations and study their long time behavior and the nature 

, of their evolutionary singularities when computation is tractable. Numerical procedures and 

5^ ' simulations are carried. 
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1 Introduction 

Structured populations are populations in which individuals differ according to variables that 
affect their reproductive or survival capacities. These variables can be phenotypic, genotypic 
or behaviorial traits that are assumed to be hereditarily transmitted from a parent to its de- 
scendants unless a mutation occurs. They can also be position, sex or age. In this article, 
we are interested in modelling the adaptive evolution of trait and age-structured populations, 
describing the ability of a population to generate and select trait diversity at an individual level 
and through the expression of the individual ages and traits. Individual-centered models allow 
us to give a realistic description of these phenomena and to obtain macroscopic equations, at a 
large population scale. 

As emphasized by Charlesworth [7], it is important to take age structure into account. Most 
populations consist in sets of individuals born over a range of past times, with behaviors that 
depend on their ages. As mentioned in [7j, there exist for many species a juvenile and an 
adult periods. The first period is devoted to maturation and growth, and the individual starts 
reproducing only during the second period. See the Tribolium population (cf. [19]) as an 
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illustration. For many species, and mammals in particular, fecundity and survival functions are 
decreasing with age (cf. [32j). Some other species exhibit long postponement of reproduction. 
The pink salmon Oncorhynchus Gorbuscha breeds at about 2 years old, which is approximatively 
its life expectation (cf. [71[22|), and the cicada Magicicada reproduces and dies at 13-17 years 
old ([7]) for instance. 

Adding an age-structure allows us to study life history traits, such as age at maturity, or to 
consider senescence phenomena, which describe the fixation of deleterious genes or the decrease 
of reproductive and survival capacities in the course of life. Many questions concerning these 
age-structured populations can be raised. How does the age influence the trait evolution? For a 
given trait, how does the probability of fixation or elimination by natural selection depend on its 
age of expression? How does the age-structure of a population infiuence the selective pressure 
exerted on individuals? Which age distributions will appear at equilibrium for a given trait? 
In this paper, our main interest is to generalize to age-structured populations the recent theory of 
Adaptive Dynamics. First models of Adaptive Dynamics have been introduced by Hofbauer and 
Sigmund ^20j, Metz et al. [26j, Dieckmann and Law [lOj and rigorous microscopic derivations 
have been obtained by Champagnat O 0], Champagnat et al. [3 [5]. In these models, the 
population is only structured by the trait of interest. 

We start with the description of a discrete stochastic individual-centered model of population 
structured by trait and age. We consider approximations of the microscopic process under a 
large population asymptotic and obtain in the limit a deterministic partial differential equation 
(PDF) involving trait and age, that generalizes classical demographic PDFs (see McKendrick 
[24J, Von Foerster [11], Gurtin and McCamy [l7|, Webb [38], Murray [28j, Charlesworth [7], 
Thieme |35j). This gives us information on the long time behavior of the population. We are 
then able to characterize rare mutation rates for which it is possible, under an additional as- 
sumption of non-coexistence in the long term of two different traits, to separate the time scale of 
ecology describing the demographic variations in the population and the time scale of evolution 
linked to the occurrence of mutations. The latter hypothesis generalizes the "Invasion implies 
Fixation" assumption of the Adaptive Dynamics theory. We prove in this case that the mi- 
croscopic process converges to the so-called Age-structured Trait Substitution Sequence Process 
that jumps from a monomorphic equilibrium (where all the individuals carry the same trait) to 
another. This process extends to populations with age-structure the Trait Substitution Sequence 
(TSS) introduced by Metz et al. [SH]. A difficulty in our case lies in the fact that the fitness 
function describing the invasion of a mutant trait in the resident population and which appears 
in the generator of the TSS process can not be computed explicitly. Taking the limit of the TSS 
when the mutation steps tend to zero gives us an age-dependent ordinary differential equation 
(ODE), generalizing the Canonical Equation proposed by Dieckmann and Law [lOj . We show 
that an explicit expression for the fitness gradient appearing in the ODE is available despite of 
the implicit definition of the fitness function for the TSS. To our knowledge, the equations that 
we establish have not been introduced in the biological literature yet and are more general than 
the ones previously proposed by Dieckmann et al. [9], Ernande et al. [E], Metz p5], Parvinen 
et al. [HU]. 

Our results point out that in the adaptive dynamics limit, the age-structured TSS and Canon- 
ical Equation involve averaged functions in age. The age-structure does not mainly affect the 
qualitative behavior of the population but plays a role in the trade-offs determining the Evolu- 
tionary Stable Strategies (ESS). With age-structure, we gain in realism and refine the results by 
including effects due to life histories and age-dependent behaviors. 

Sections [2] to [3] are devoted to the presentation of the three models of interest: the microscopic 
model, the TSS model, and the Canonical Equation. Then, we state the limit theorems linking 
these models (Theorems 13.11 and 14. ip . Several examples are considered, in Section [5l we begin 
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with a logistic model with senescence in the birth rate and size dependence in the competition 
term. The long time behaviors of the TSS process and of the solution of the Canonical Equation 
are studied. We highlight differences with the case without age-structure. In our model, senes- 
cence acts as a penalization term that favors fast reproduction to growth. We then consider a 
model belonging to a class that we call age-logistic. Computation becomes rapidly complicated. 
In Section [6] we investigate an example where the competition kernel is nonlocal, asymmetric in 
trait (cf. Kisdi [23j) and decreasing in the age of the competitors. The explicit expression of the 
fitness gradient allows us to compute the ESS value and to show that coexistence is possible. 
The existence of an "evolutionary branching" is suggested by numerical simulations. 

2 Microscopic model 

We start with a stochastic microscopic model of age and trait structured asexual population 
describing the dynamics at the individual level. We take into account births, either clonal 
or with mutation and death, natural or due to the competition with other individuals. A 
large population limit is studied. We provide by this way a microscopic justification of PDEs 
generalizing models introduced in Demography ([71 dU ESI EH])- This gives us a better 
understanding of the large time behavior of the microscopic process and allows us to define the 
rare mutations asymptotics leading to Adaptive Dynamics limits. To our knowledge, there are 
only a few models dealing with trait and age-structured populations, all deterministic (Rotenberg 
[M| , Mischler et al. [27] ) and the stochastic models taking age into account do not consider trait 
evolution, see Kendah [21], Athreya and Ney [l], Doney [12|, Oelschlager [29], Tran [37]. We 
generalize all these models by introducing a dependence between trait and age at the individual 
level and by taking into account mutation and competition between individuals, which yields 
nonlinearity in the limiting phenomena. 

2.1 Microscopic description 

Individuals are characterized by a trait x belonging to a compact set X of and by their 
physical age a G M+. We set X := X x The population is discrete and described by the 
point measure 

{Zul) 

Zt{dx,da) = J2 \x,{t)Mt))- (^-l) 
1=1 

Each individual is represented by a Dirac mass on its trait and age, and {Zt, 1) is the population 
size at time t. We denote by Aipi^) tbe set of finite measures on X. For Z G Mf{X) and for 
a real- valued bounded measurable function /, we define {Z, f) = fdZ. 

An individual with trait x G X and age a G M+ in the population described by Z gives birth to 
a new individual with rate b{x, a). With probability p G [0, 1] the new individual is a mutant of 
age and trait x + h, where h is drawn from a probability distribution k{x, h) dh with support 
on X — {x} = {y — x\y £ X} (so that x + h £ X). With probability 1 — p £ [0, 1], the new 
individual is clonal, with age and trait x. The individual dies with rate d{x, a) + ZU{x, a) 
where d{x,a) is the natural death rate and ZU{x,a) = f^U{{x,a),{y,a))Z{dy,da). Here, U 
is a real-valued kernel describing the interaction exerted by an individual with trait y and age 
a on an individual with trait x and age a. During their life, individuals age with velocity 1, so 
that the age at time t £ R_|_ of an individual born at time c < t is a = t — c. (Notice that the 
birth date c varies from an individual to another). 

Let us now describe the generator of the A4ir(A')-valued Markov process Z. This generator L 
sums the aging phenomenon and the ecological dynamics of the population. Let C^'^(A',M) be 
the space of continuous bounded real- valued functions on X with bounded continuous derivatives 
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with respect to the age variable. As developed in Dawson [H] Theorem 3.2.6, the set of cylindrical 
functions defined for each ^^ G Mf{^) by Ff{fi) = F{{n, /)),^with F G Cl{R) and / G C°'^(^), 
generates the set of bounded measurable functions on A4f{X). For such function, 

LFfifi) = {f,,dafi.))F'fif,) 

+ £ [{d{x, a) + nUix, a)){Ff{fi - 6^^^a)) - Ff{ii)) + a)(l - + 5(^,o)) " ^/(m)) 

fi{dx,da). (2.2) 



+ 6(x, a)p / (Ffifi + 5(^+/j,o)) - Ff{n))k{x, h)dh 
In the sequel, we make the following assumption. 

Assumption 2.1. The functions b, d and U are assumed to be of class and x ^ k{x,h) 
is Lipschitz continuous, uniformly in h. We suppose that there exist strictly positive constants 
b, d, d, U, U such that V(x, a) G X, V(y, a) G X , 



U{{x,a),{y,a)) < U, / {sup k{x, h))dh < 



0<b{x,a)<b, d<d{x,a)<d, U_< U{{x,a),{y,a)) < U, / {sup k{x, h))dh < +oo. 



Under Assumption 12. H the process Z can be obtained as the unique strong solution of a 
stochastic differential equation driven by a multivariate Poisson point measure corresponding 
to the dynamics described above (we can adapt Fournier Meleard [15] and Tran [37] concerning 
respectively trait-structured and age-structured cases). 

Let us mention different forms for the birth rate that can be found in the literature. The 
period devoted to reproduction is very different from one specie to another. For some species, 
there is no reproduction in the first growth period (as juvenile or larva states, cf. Henson |19]). 
The parameters of interest are the maturity age cm and the constant birth rate bi, which are 
the heritable traits. The birth rate for an individual with traits {aM,bi) and age a can be 
b{aM,bi,a) = bil^^y^^^^y 

Many other species reproduce early in life, and the reproduction function is decreasing with 
age. The heritable traits submitted to mutation are then the initial reproduction rate bi and 
the senescence parameter 62. As in (cf. Webb [38] p. 39 Ex. 2.1), the birth rate of an individual 
with traits (61,62) and age a can be 6(61,62,0) = bie~^'^"-. 

Conversely, for other species ( |38j p. 41 Ex. 2.2), reproduction happens at the end of life and 
the birth rate of an individual with age o can be modelled as 6(61,62,0) = 6i(l — e"''^"). 
The individual trait will be (61,62), with 61 the reproduction rate at maturation and 62 the 
maturation parameter. 

The competition kernel U describes the competition intensity between two individuals and can 
depend on their respective ages and traits. The simplest case is the (density dependent) logistic 
one, for which ?7((x, a), (y, a)) = r]{x,a). The interaction exerted on an individual is then 
proportional to the total number of individuals in the population, and thus nonlocal. Another 
interesting case is when U{{x,a),{y,a)) = r]{x,a)U{x,y). Then interaction is nonlocal in age 
and local in the trait space. Such population will be called age-logistic. 

In Sections 5 and 6, we develop examples of asexual age and size-structured populations. We first 
consider a logistic population with senescence in the birth rate. Computation can be developed 
in details. We compare it with the corresponding model without senescence and with a model 
with age-logistic interaction. In the last model, we investigate a case for which the competition 
kernel is local in trait and age. 
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In these examples, the trait space is the size space [0, 4] and the size trait is heritable except 
when a mutation occurs. With probability p g]0, 1[, the descendant of an individual with trait 
a; is a mutant with trait x + h. Mutation amplitude h is distributed on [—x, 4 — x] following the 
centered law k{x,h)dh, so that x + h remains in [0,4]. With probability 1 — p, the descendant 
is clonal. 

Example 1: The individuals reproduce at rate 

6(x,a) =x(4-x)e-°. (2.3) 

Dependence in x explains as follows: if the individual is small, it has not enough energy to 
reproduce with high rate. If it is large, descendants' size will also be large by hereditary trans- 
mission, and more energy is required to produce such descendants than smaller ones. The term 
e~" expresses a senescence phenomenon and the reproduction rate is decreasing with age. The 
logistic death rate of an individual with trait x is given whatever its age by 

d{x, {Z, 1)) = \ + 0.001 • (4 - x){Z, 1). (2.4) 

The term 1/4 is the natural death rate and the logistic term has intensity 0.001(4 — x), meaning 
that big individuals are less sensitive to competition. The trait maximizing the birth rate is 
X = 2. Since the logistic competition favors individuals with size close to 4, we can expect that 
the optimal traits, taking birth and death rates into account, are lying in the interval [2,4]. 
In Section 5, computation will allow us to obtain the value of the trait equilibrium (evolutionary 
stable strategy) in which the population will stabilize in the evolution time scale. We will 
compare this model with the similar one with b{x) = x(4 — x), (a size-structured population 
without age structure) and will show that the senescence phenomenon appears as a penalization 
factor, reducing the equilibrium trait value. We also compare the first model with the one in 
which the logistic interaction term 0.001(4 — x){Z, 1) is replaced by 

a [ U{x,y)Z{dy,da), where U{x,y)=c(l-— (2.5) 

Jx \ l + uexp{-k[x -y))J 

is Kisdi's asymmetric competition function [23j (C = 0.002, z/ = 1.2 and k = 4). Its sigmoid 
shape models the fact that the competition on an individual is mainly exerted by larger ones. 
We see how intricate computation can be. 

Example 2: We investigate a model where the competition kernel is a function of the traits and 
ages of the competitors. An individual with trait x and age a reproduces at rate b{x) = x(4 — x). 
There is no natural death and the action exerted by a competitor of trait y and age a is de- 
scribed by the interaction function rj(a) (1 + e~") U{x,y) of separable multiplicative form, with 
U defined in (j2.5p . The term 1 + e~" is a senescence parameter diminishing the interaction 
intensity with the competitor's age. When is a constant function, we recover models studied 
by [21 El] with a deterministic approach. In Section 6, we investigate the case r/(a) = a. 

Simulations using an individual-based algorithm are given in Figure[I](c). They present the trait- 
support of the point measure Z defined in (|2.ip for the three examples previously described. We 
remark a qualitative difference between cases (a) and (b) and case (c). In the two first plots, the 
trait-support is very dense. Conversely, in Figure [1] (c) , the population separates into smaller 
groups once it has reached a certain trait value which seems to be approximatively 3.2. This 
phenomenon is known as branching in the biological literature. 
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Figure 1: Simulation of the microscopic process, using an individual- based algorithm: the traits 
in the population are represented in function of time. We start with n = 2000 individuals. The 
initial traits are drawn uniformly between and 1.3, and the ages are drawn in an exponential 
distribution of parameter 2. For the mutations, p = 0.13 and k is the gaussian law with mean 
and variance 0.15, conditioned to stay in [0,4]. (a): Logistic population with age structure 112.31 
\2.4^ . (b): Logistic population without age structure and \2.4^ . (c): Age- structured population 
with nonlocal interaction kernel (Example 2). 



2.2 Large population and rare mutation renormalizations 

We are interested in studying approximations of the age-trait dynamics presented in Section 12.11 
under a large initial population and rare mutation assumption. 

The processes of interest are the renormalized population processes given by: 

yneN*,yt€R+, Z^{dx,da) = - ^ <J(., (2.6) 

1=1 

The assumptions are as follows. 

Assumption 2.2. (i) 3e > 0, sup„gi^, E {{Z^, 1)2+^) < +oo. 

(a) The sequence Zq converges in law in Aipi'^) to a finite measure ^o- 

(Hi) The ecological parameters stay unchanged, except the interaction kernel Un = This 
hypothesis corresponds heuristically to the expression of a resource constraint. The death rate of 
an individual (x,a) of Z^^ is thus equal to d{x,a) + Z^U{x,a). 

(iv) The mutation probability is given by UnP, with Un, p £ [0, 1] and Un is assumed to decrease 
to when n +00. 

Adapting tightness-compactness results in Fournier and Meleard [T5], one can prove that the 
sequence converges to a deterministic limit when n -\-oo. More precisely: 

Proposition 2.3. (cf. Tran ^36] (Section 3.2)) Under Assumptions \2.1\ and \2.^ the sequence 
converges in D(]R+, Al^T^f A")) to a deterministic continuous process ^ G C{M.^,Mf{X)) 
characterized as the unique solution of the evolution equation: V/ G C^'^ {X ,R) , \/t G 

iCtJ) = {Co,f)+ [ [[daf{x,a) + f{x,Q)b{x,a) - f{x,a)(d{x,a) CsU{x,a))]^s{dx,da) ds, 

JQ JX 

(2.7) 
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Let us roughly give the main ideas of the proof. Since every individual is weighted by the 
size of population jumps converges to zero, and the limiting values of are continuous 

processes. Moreover, for / G C^'^(A',M), {ZJ^,f) writes as the sum of a finite variation term 
and a random term whose square is of order 1/n. The random part disappears when n — > +00 
and the limit is deterministic. In the finite variation part, mutations also disappear in the limit, 
since (un)neN* tend to 0, leading to (12. 7p . No diversity appears since there is no mutation. 
Uniqueness of the solution of (12. 7p implies uniqueness of the limiting value of (Z")„gN*. 

In cases where the support of the trait-marginal of ^0 is a singleton {xq}, respectively a pair 
{xo,y}, the population remains monomorphic, respectively dimorphic. Equation (j2.7p is then 
parameterized by these values and relations between ()2.7p and classical partial differential equa- 
tions have been studied by Tran ([37] Propositions 3.4 and 3.6). 

Proposition 2.4. (i) Assume that S,o{dx,da) = Sxo{dx)mo{xo, a)da. Then for every t > 0, 
£,t{dx,da) = 6xo{dx)m{xo,a,t)da and the function m{xo,a,t) is the unique weak function 
solution of the partial differential equation parameterized by xq : 

dm dm if \ 

-^(xo, a, t)+ — {xQ, a,t) = - \^d{xo, a) + J U ((xq, a), (xq, a))m(xo, a, t)daj m(xo, a, t), 

r+00 

m(xo,0, t)= / b{xo,a)m{xo,a,t)da, m(xo, a, 0) = mo(xo, a). (2-8) 



(a) Assume that (,o{dx, da) = 6xQ{dx)mQ{xo, a)da + 6y{dx)mo{y, a)da. Then for every t > 0, 

(,t{dx, da) = 6xo{dx)m{xo, a, t)da + 6y{dx)m{y, a, t)da 
where {m{xo,a,t),m{y,a,t)) is the unique weak solution of 

dm dm 

— [xo,a,t) + —[xo,a,t) = 

d{xo, a) + I {U{{xo, a), (xq, a))m(xo, a, t) + U{{xq, a), (y, a))m{y, a, t)) da ] m{xo, a, t) 



dm , , dm , , 
^(y,a,t) + ^(y,a,t) = 

a) + / (^((y. o)) (xo, a))?n(xo, a, t) + C/ ((y, a), (y, a))m{y, a, t)) da ] m{y, a, t), 
Jr+ J 

f+oo r+oo 

m(xo,0,t)= / b{xo,a)m{xo,a,t)da, m{y,0,t) = / b{y,a)m{y,a,t)da, 
Jo Jo 

m(xo,a,0) = mo(xo,a), m{y,a,0) = mo{y,a). (2.9) 

Equation (12. 8p generalizes Demography equations, as McKendrick-Von Foerster or Gurtin- 
McCamy Equations (see [MlOllTi). 

In the sequel, we will make the following assumptions concerning the long time behavior of the 
solutions of I^M and UiTM. 



Assumption 2.5. Let xq £ X and mQ{xQ,a)da G A^i?(M+). The solution of h2. admits a 
unique nontrivial stable stationary solution iri{xQ,a) such that m{xo, a,t)da converges for the 
weak convergence topology in M.fO^+) to in{xQ, a)da when t — > +00. We will denote by 

r+co 

Mxo = / m{xo,a)da (2.10) 



the mass of the stationary age measure and by the space-time stationary measure S^xoidx, da) 
5xi^{dx)m{xQ, a)da. 
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Under Assumption 12. 5[ m(xo,a) is solution of 

din ^ ^ 

— {xo,a) = -d{xo,a,xo)m{xo,a), d(xo,a,xo) ■.= d{xo,a)+ I U{{xQ,a),{xo,a))m{xo,a)da, 
oa ' 



R+ 



(2.11) 



with the boundary condition: 



r+oo 

fh{xQ,0) = / b{xQ,a)rh{xo,a)da. (2-12) 
Jo 

Solutions to ()2.1ip have the form 

■m{xo,a) = m{xo,0) exp ^— J d{xo,a,XQ)da^ . (2.13) 
From (|2.13|) and (j2.12p . we obtain the following balance condition 

b{xo, a)e- i'o °° 'i(^o,",^o)da^^ ^ ^2.14) 

A necessary condition to get Assumption 12.51 is then that 

r+oo 

i?o(xo):=/ 6(xo,a)e-^o''««{^o,«)rf"^Q > 1_ (2.15) 
Jo 

The term Rq{x()), called net reproduction rate, is the integral in age of the birth rate weighted 
by the survival probability in absence of competition. It's the well known threshold between sub 
and super-criticality in age-structured models (see [H 111 [TTl [38] ) . If the birth and death rates 
b{x) and d{x) do not depend on age, Ro{x) = b{x)/d{x) > 1 if and only if b{x) > d{x), which is 
the standard super-criticality condition. 

With the notation of Section 12.11 in the age- logistic case and for a monomorphic population with 
trait xq, the death rate in (12. 8p equals d{xo,a) + r]{xo,a)U{xo,xo) J^^ m{xo,a,t)da. Equation 
(12. Sp is then a Gurtin-McCamy equation parameterized by xq and Condition (I2.15P is also 
sufficient to obtain Assumption 12.51 as proved in Webb [38], Section 5.4. 

Let us now introduce assumptions describing the "invasion implies fixation" principle, for the 
age-dependent system defined by ()2.9p . Assertion (ii) impedes the co-existence of the two traits 
Xq and y in the long time. 

Assumption 2.6. (i) Let Xo,y S ^ with Ro{xq) > 1. One of the following assumptions holds: 

Either: [ 6(xo, a)e- i'o^(^»'"'2')'^°da > 1 and [ 6(1/, a)e" i'o^^^'^'^o^'^^da < 1 (2.16) 
Jr+ Jr+ 

Or: [ b{xo,a)e-^Sdi.^^^^,y)d<^da<l and I 6(y, a)e- ^(^'"'^"^'^"da > 1, (2.17) 



where 



-co 



d(x, o, y) := (i(x, o) + / U{{x,a),{y,a))7n{y,a)da (2.18) 







is the death rate of an individual {x,a) in the population at equilibrium ^y. 

(ii) The solution of \2.1^ with the dimorphic initial condition ^^[dx^da) = 6xQ{dx)mQ{xQ,a)da + 
5y{dx)mo{y,a)da converges as t ^ 00 to if (KIW ^•s satisfied and to if 1^2.17^ is satisfied 
(for the weak convergence topology in Mf{X)). 
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For a logistic dimorphic population with traits xq and y {Rq{xo) > 1), Assumption [JS] is fulfilled 
as soon as U (xq, xo)U{y, y) — U{xo,y)U {y, xq) < and one of the following assertion is satisfied: 

U{xo,xo)M^,-U{xo,y)My <0 and U{y,y)My - U{y,xo)M^^ > 0, 
oiU{xo,xo)M^,-U{xo,y)My>0 and Uiy,y)My - U{y,xo)M^o < 0. 

(Let us recall that M^q has been defined in (I2.10p ). 



3 Age Structured Trait Substitution Sequence Process 

In order to obtain the Adaptive Dynamics limit, we consider here the same mutation scale as in 
the work of Champagnat [3] without age-structure. We assume that 

V y > 0, exp(—Vn) = oiun) and Un = o ( — | . (3.1) 

\ n log n / 

This scaling is here again the right one to derive from the microscopic process a jump process 
generalizing the Trait Substitution Sequence (see Metz et al. [26], Champagnat [4J) for age- 
structured populations. It is obtained from a fine study of the different time steps in the invasion 
process of the resident population by a mutant one. This study has been done by Champagnat 
in [1] and adapted to include age structure in Tran's thesis |36) . Let us roughly recall the main 
ideas justifying (jS.ip . They are summarized in Figure [2j If no mutation occurs, then with a 
probability that tends to 1 when n — > -|-oo, a monomorphic population with initial size n enters 
a given neighborhood of its equilibrium after a sufficiently large time that does not depend on 
n. By large deviation results, the microscopic process stays in this neighborhood during an 
exponential time e"^. Mutation time is of order \/{nUn), and an invasion period is proved to 
be of order logn. Thus the condition Un = o(l/(nlogn)) in (j3.ip implies that the mutations 
are sufficiently rare so that the population has returned to its monomorphic equilibrium when 
the next mutation occurs. The condition exp(— T/n) = o{un) tells us that the mutations occur 
sufficiently often so that a mutant appears before a rare event drives the resident population far 
from its equilibrium. Hence, if we change time and consider (Z"j.^^ pngN* at the mutation time 
scale, we will obtain in the limit, when n +oo, a process where the transition periods have 
disappeared and where only the sequence of equilibrium states remains. Under our hypotheses 
preventing the trait co-existence, the latter have a singleton trait-support and then reduce to 
age-measures. Therefore, the limiting process jumps from an age-measure to another, each of 
these parameterized by a trait. This provides a generalization of the TSS proposed by Metz et 
al. [26]. The result is stated in Theorem 13.11 and proved in Appendix A. 

Theorem 3.1. Let us consider [Z"-) defined by \2. 6\) with a monomorphic trait support {Xq} 
and satisfying Assumptions \2.1l \2.2\ with S,o{dx,da) = 6xoidx)mo{XQ, a)da. Assume moreover 
that {un)nen satisfies 113. 1\) and that Assumptions \2.5\ and \2.6\ are satisfied for each xo,y G X 
with Rq(xo) > 1. Then, for every t S M.^, the sequence (Z^^^^ pneN* converges in law in 

Aipi^) to the measure Zt defined by 

Zt{dx,da) = Sxt{dx)fh{Xt,a) da, (3.2) 
where X is a Markov jump process of generator defined for (f) G Bb{X, M), x ^ X by 



L(t){x) 



{(j){x + h) — (j)[x)) p [ I h{x,a)m{x,a)da\ {1 — zq{x + h,x)) 



k{x,h)dh, (3.3) 
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population size 




Figure 2: Large time behaviour of the microscopic process Z" and of its deterministic approximation 
^. On compact time intervals and for sufficiently large n, the behaviour of follows the one of its 
approximation ^. On a large time scale, if no mutation occurs, Z" leaves the neighborhood of the sta- 
tionary equilibrium of ^ at T" after an exponentially long period and drives the population to extinction. 
However, before it gets extinct, the population can be invaded by a mutant in a time scale of order 1 /nUn ■ 
Under the "Invasion implies fixation" Assumption \2. b\ if the mutant population does not die, it replaces 
the resident population after a transition period of order logn. 

where ZQ{y,x) is the smallest solution in [0,1] of: 

f+oo ^ 

z= / e^^-^^/o" a, x)e-/o"'^(2/'«'^')'i«c^a. (3.4) 

Jo 

The death rate d{y, a, x) has been defined in i2.18\) . 

Corollary 3.2. The sequence (Zy^^^^^^^p^* converges to Z in B(M+, in the sense of 
finite- dimensional distributions. 

The process X is called Age Structured Trait Substitution Sequence Process. 

We have already given the main ideas of Theorem I3.1f s proof. Let us comment on the different 
terms of the generator ()3.3p . 

1. The term pk{x, h) b{x, a)fh{x, a)da is the rate at which the monomorphic population of 
trait X, in its equilibrium, generates a mutant of trait x -\- h. This term equals pb{x)Mxk{x, h) 
in the absence of age-structure. 

2. The term (1 — zo{x-\- h, x)) is called fitness of the mutant trait x-\-hm the population of trait 
x and is the probability that the mutant descendants invade the resident population. In most 
cases, this probability zq{x + h, x) can not be computed explicitly conversely to the case without 
age-structure. Nevertheless, its implicit definition given in Equation ()3.4p makes possible its 
numerical computation as developed in the examples of Sections [5] and El 

Let us now focus on the way we establish the equation defining the fitness function ZQ{y,xo) 
where xq £ is the resident trait and y the mutant trait. When the mutant y appears in 
the monomorphic resident population with trait xq, we can neglect in a first approximation 
the mutant population and the deviations of the resident population from its equilibrium. The 
mutant population is then compared to a linear age-structured birth and death process with 
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parameters b{y,a) and d{y,a,XQ). For this process, the extinction probabihty is the smallest 
solution ZQ^y^xo) in [0, 1] of ()3.4p . Next proposition proved in Appendix lAl identifies zo{y,xo) as 
the extinction probability of the real mutant progeny in large population. 

Proposition 3.3. Consider the process as in Theorem \3.1\ and starting from Z'Q{dx,da) = 
6xo{dx)qQ{xo,da) + ^6(^y Qj{dx,da), where 6xo{dx)qQ{xo,da) is a point measure weighted by 1/n 

with support in {xq} x M+ and converging to S^xq in ■Mf{'^)- We denote by ^xoq^y ^^"^ 
Let T be the first time at which a mutation occurs, 6 the first time of return to a monomorphic 
population and V the survival trait at time 9. Then, under the assumptions of Theorem \3. 11 



Figure 3: An age- structured process in continuous time and its discrete time generation-tree. 

A proof of ()3.4p is adapted from Doney [12]. Noticing that the extinction of the mutant progeny 
for the linear age-structured birth and death process is equivalent to the extinction of the 
discrete-time Galton- Watson process corresponding to the underlying generation process (see 
Figure[3]), zo{y, xq) is characterized as the smallest fixed point on [0, 1] of the generating function 
of the number of children and satisfies: 




time 





(3.5) 



(3.6) 



Equation (13.60 is obtained from ()3.5p by integration by parts. 



A main issue is to know whether extinction happens almost surely or not. 




(3.7) 



then zo(?/,Xo) < 1 and defines a function in both variables. 

(a) Else, zo{y,xo) = 1. 

(Hi) Finally, when y = xq, zo{xq,xo) = 1. 
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Of course, z = 1 is an obvious solution of (j3.6p . Each other solution is a zero of F where 

f + OO ^ 

Fiz, y,x)= / b{y, a)e^^-^^ Ky,^)du-Jo d{y,u,xo)du^^ _ -^^ ^3 



The proof is therefore an immediate consequence of the following lemma. 

Lemma 3.5. For each pair of traits x,y, the equation F{z,y,x) = admits a unique solution 
z = g{y,x) on R+ of class in both variables and such that g{x,x) = 1. 

Proof. The function F is of class with a positive partial derivative with respect to z. It is 
thus strictly increasing in z. Since ff^°° b{y,a)e~ ■f^o Ky^u)du^^ _ ^- b{y,u)du _ i ^ then 
F{0,y,x) < 0. Moreover liniz^+oo F{z,y,x) = +00. We deduce that for every x,y £ X, there 
is a unique solution z = g{y, x) on M+ to F{z, y,x) = which is equal to 1 if and only if 

fo~°° ^^y^ ^)^~ '^'^^'"'^"•''^"(ia = 1. This condition is true for y = xq by ()2.14p . Using the Imphcit 
Function Theorem we obtain moreover that g is of class C^. ■ 

In both cases, we can notice that: 

zo{y, xq) = g{y, xq) A 1 and 1 - zo{y, xq) = [1 - g{y, xo)]+ , (3.9) 
where [.]+ denotes the positive part. 



Remark 3.6. For a population without age-structure Equation ()3.5I) can be solved explicitly 
and we recover the infinitesimal generator introduced by Metz et al. [26j with: 



1 - zo{y,xo) 



b{y)-d{y)-Uiy,xo)M,, 



(3.10) 



4 Canonical Equation for an Age-structured population 

We are interested in the limit of the Age-structured Trait Substitution Sequence Process when 
the mutation step tends to zero. Generalizing the approach of Dieckmann and Law [10] and 
Champagnat [3j in the case without age, we consider, for e > 0, the infinitesimal generator 
defined for every (/> e i3f,(A',M) and x £ X by: 



F 



'(f){x) = ^ [ {(j){x + eh) — (j){x)) p [ b{x,a)m{x,a)da{l — zq{x + eh,x))k{x,h)dh. (4.1) 



When £ ^ 0, the sequence of such renormalized TSS-processes converges to the solution of an 
ODE that generalizes the Canonical Equation introduced by Diekmann and Law [lOj: 

Theorem 4.1. Under Assumptions \2. 1[ \2. 51 \2. 61 the sequence (X^)£>o converges in probability, 
for the Skorohod topology on B (IR+, A") to the solution of the following ODE: 

-J- = P j b{x,a)m.{x,a)da f hD\zQ{x,x)k{x,h)dh, (4-2) 



dx \x{t + h) - x{t)\ , / X 1- zo{x,x) - zo{x + £h,x) 

where —— = hm ; , ana DhZnix.x) := lim 

dt h^o h £ ^ ^ 

e > 
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The Proof of Theorem I4.H based on a tightness-uniqueness argument, can be adapted from 
Theorem 1 in Champagnat [3]. 

Let us remark that D\zq{x,x) is nonnegative since zq{x + eh,x) G [0,1] and zo{x,x) = 1 by 
Proposition I3.4[ Using (|3.9|) . one gets 



Dlzo{x,x) 



Mm 

£ ^ 
£ > 



1 — g{x + eh, x) 



(4.3) 



Hence, as in the classical case, the evolution follows the directions where the fitness gradient 
Dj^zo{x,x) is positive, and along which the extinction probability y >—>■ zo{y,x) is a decreasing 
function. 

Although zo is implicitly defined, an explicit expression of Dl^zo{x,x) can be established. 

Proposition 4.2. Let us consider the scalar case d = 1 for sake of simplicity. Under ^s- 
sumptions \2.1\ [275[ \2.6[ for h > 0, D'Izo{x,x) = h[—dig{x,x)]- and for h < 0, D'Izq{x,x) = 
h[—dig{x,x)]-^-, with g defined in Lemma VJ. 51 and 



dig{x,x) 



dib{x, a) — b{x, a) 



a 


+ 00 



did{x,a) + / diU{{x,a),{x,u))fh{x,u)du\ da 



X e-io"'^(^'".^)'i°cia X ( ^ b{x,a)i^J b{x,a)da] e' '^i^''''^)'i"da] . (4.4) 



-1 



Proof. We replace y by x + e in ()3.8p and consider the expansion with respect to e, using: 

b{x + £,a) = b{x, a) + £dib{x, a) + o(e), 
d{x + e, a) = d(x, a) + edid{x, a) + o(e), 

U{{x + £,a), {x,a)) = U ((x, a), {x, a)) + £diU{{x, a), (x, a)) + o(e), 
(^(x + e, x) = 1 + £dig{x, x) + o(e). 

Identifying the terms of order £, we obtain (14. 4p . ■ 

Remark 4.3. (i) When c?ig(x,x) > 0, DIzq{x,x) = and L'i]^2:o(3;) 3;) = dig{x,x). When 
dig{x,x) < 0, DIzo{x,x) = —dig{x,x) and Dl^iZo{x,x) = 0. 

(ii) In the scalar case without age-structure, ()4.4|) allows us to recover the expression of the 
classical fitness gradient: 



dIzo{x,x) 



b'{x){d{x)+U{x,x)M^)-bix){did{x) + d2U{x,x)M^) 

6(^)2 

^b{y)-{d{y) + U{y,x)M^) 



Ky) 



y=x 



5 Example 1 

We now present several examples for which we specify the adaptive dynamics approximations 
and study their behavior. We used R for simulations which illustrate our purpose, and Maple 
for formal calculus, when computation becomes too technical. 
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5.1 A logistic age and size-structured population 



We develop in this section a simple example highlighting the difficulties that appear when 
considering a trait and age-structured population and for which computations can be carried 
explicitly. The birth and death rates that are used have been specified and explained in (12. 3p 
and (|2.4|) . Compared with models with only trait structure, the specificity here lies in the 
introduction of a senescence term. The latter, even if it is simple, introduces a notion of life 
history that will have an effect on the traits selected through evolution. 

5.1.1 Monomorphic equilibrium in large populations 

The large population approximation is given by Proposition 12.41 



dm , , dm 
— {x,a,t)+—{x,a,t) 



- + 0.001(4 — rc) / m{x,a,t)dajm{x,a,t) 
m{x,0,t) = I x{4: — x)e~"'m{x,a,t)da. (5.1) 





From ()2.15p . there exists a non trivial stationary solution if and only if: 

" + 00 



Ro{x) = I ~ x(4 - x)e-5'^/^da > 1 ^ ^""^^ > 1 
Jo 5 

|0.35,3.65[. (5.2) 



x e 



2-^,2 + ^ 



Under ()5.2p . Assumption 12.51 is satisfied. 

Proposition 5.1. // < [5.^|) is satisfied, then h5.1\) admits a unique nontrivial stationary solution 



(a;(4 - x) - 5/4) (2;(4 - x) - 1) , , , ^ x s , 
m{x, a) = omi{4 - x) ("^^^^ " " ^^^^ " ^^'^^ 



Proof. Any stationary solution is of the form m(x, a) = m(x, 0) exp ( —(1/4 + 0.001(4 — x)Mx)a 
by ()2.13p . Plugging this expression into ()2.12p gives 



m(x, 0) =m(x, 0) /^^ x(4 - x)e-(t+0-00i(4-x)M.)a^^ ^ ^^^^ 0^ ^i^^ ^) _ _ 4) 

Jo I + 0.001(4 - x)Mx 

Since we are looking for a nontrivial solution, m(x, 0) 7^ and necessarily: 

0.001(4 -x) ■ ^ ■ ^ 

The definition of implies = m(x, 0)/(l/4 + 0.001(4 - x)M^) which gives m(x, 0). ■ 

5.1.2 Invasibility 

With our choice of competition kernel. Assumptions 12 . 6l are satisfied. Let us consider the invasion 
phenomena. An explicit expression can not be obtained for ZQ{y,x), but we can compute it 
numerically (see Figure S]), which allows us to simulate the paths of the TSS process (see Figure 
[S]). Even if zo(y,x) remains implicit, the study of the domain of invasibility, that consists in 
pairs (x, y) of traits y with positive fitness in a resident population of trait x, can be carried 
explicitly. It brings information on the long time behavior of the Age structured TSS process. 
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Proposition 5.2. Let x G [0,4] satisfy 15. The mutant traits y which can invade the 
monomorphic resident population with trait x belong to ] min(x, /(x)), max(x, where: 

/(x):=4--^. (5.6) 

4 — 2; 

Proof. From Proposition 13.4^ we know that ZQ{y,x) < 1 (i.e. the invasion of the resident 
population by the mutant one is possible) if and only if /q^°° b{y, a)e"'^(^''*^^^"(ia > 1. Since 

»+oo — r+oo — — — 

b{y,a)e-'^^y'^^^'^''da= / a)e-'^(2''*^«)''+('^(^'*^«)-'^(^'*^-)''da, (5.7) 



and 

r+oo 

b{y,a)e-'^^y'^^yUa = l 



by the balance condition (I2.14p . Equation ()5.7p is satisfied if and only if d{y, My) — d{y, M^) > 0. 
This is equivalent to: 

M, - M. > ^^ii^^l^ < 

^ 0.001(4 -x) 0.001(4 -y) 

5 5 
4^x(4 - x)(4 - y) - -(4 -y) < y(4 - x)(4 - y) - -(4 - x) 

5 

4^{x - y)(4 - x)(4 -y)--(4-y-4 + x) <0 



^(x - y) 
Case 1: if x > y (j5.8p becomes: 



(4-x)(4-y)-^ 



< 0. (5.8) 



rA N 5 , 5/4 

(4 - x)(4 - y) < - <;4>y > 4 



4 (4 - x) ' 

and hence /(x) < y < x, with / defined in (|5.6p . 

Case 2: if x < y, we obtain with similar computation that x < y < /(x). ■ 

Corollary 5.3. Let xq G [0,4] 6e an initial condition that satisfies \5.S^) . The TSS process 
{Xt)t&^^ starting from xq converges almost surely to the unique fixed point x* of f on [0,4] 
when t +00. 

Proof. The function / is a decreasing concave function satisfying 

Vx G [0, 4[, / o /(x) =4 7^^^^ = 4 - = ^- 



i-x ) 4-x 

Let's remark that / has a unique fixed point in [0,4] given by 

X* :=4- ^ « 2.88. (5.9) 

If X < y < /(x), then x = / o /(x) < /(y) < /(x), and if x > y > /(x) then x = / o 
/(x) > /(y) > /(x). Thus, the new interval of invasible traits [min(y, /(y)), max(y, /(y))] C 
[min(x, /(x)), max(x, /(x))]. 

Hence the sequence [min(Xt, /(X^)), max(Xt, /(X^))] is almost surely a strictly decreasing se- 
quence for the inclusion and it converges to x*. Since a monomorphic population with trait 
X* cannot be invaded, x* is an Evolutionary Stable Strategy (ESS) in the sense developed by 
Diekmann [11], Metz et al. [26j and Geritz et al. [16]. ■ 
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(a) (b) 




Figure 4: (a): Probability of extinction as a function of the resident and mutant traits. For every 
possible values of resident and mutant traits on a grid in [0,4], we have computed the solution 
of 113. 5\) numerically, (b): Graph of function f defined in \5.b]) . It gives the regions of the left 
figure where ZQ{y,x) = 1 and zo(y, x) < 1. This graph is called Pairwise Invasibility Plot (PIP). 

5.1.3 Simulations of the Age-structured TSS 

The Age- structured TSS Z defined by ()3.2|) is a measure-valued process jumping from an equi- 
librium measure to another. Each equilibrium measure is characterized by a trait x and by the 
density 7n{x,a) defined in ()5.3p . Simulations are given in Figure [5] below. 

(a) (b) 




Figure 5: Age- structured Trait Substitution Sequence process: the measure-valued process jumps 
from one equilibrium age-distribution to another one. (a): Each strip describes an equilibrium 
m{x,a) as age-function. These equilibria are given in ( [5.3|) . (h): the trait-valued process X 
describing the trait jumps. 

On the Figure [5l the initial resident trait is xq = 0.552 and the size of the population at its 
equilibrium is M^^ ~ 189.47. In the beginning, successful invasions favours increasing traits and 
Mx increases too, as well as the jump rate. For instance, after 4 successful invasions, the trait 
is approximatively x = 2 and M2 = 1375. When the trait approaches x* ~ 2.88, the interval of 
possible invading traits decreases and the extinction probability of possible successful mutants 
tends to one making invasions rare. 
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5.1.4 Age-structured Canonical Equation 

When the mutation step decreases to zero and time accelerates as in Section [H the evolution is 
described by the age-structured Canonical Equation (j4.2|) . 

Proposition 5.4. The fitness gradient DIzq{x,x) is obtained from ( |^.^ with: 

(4x2 _ 322; + 59)(2;2 



and the Canonical Equation has the following explicit form: 
dx 



(5.10) 



dt 



p[ [dig{x,x)]_ / h k{x,h)dh — [dig{x, x)], / h k{x,h)dh 



{-x^ + 4x - 5/4)(-x2 + 4x - 1) 



.ip(x). (5.11) 
0.001(4 -x) v>v ; V ; 

Proof. The fitness zq{x + e, x) is the smallest solution in [0, 1] of Equation (13. 6|) . which writes 

(5.12) 



z-l={z-l) (x + e)(4-x-e) 
Jo 



For the case where 

f + OO 



1 > 



5(x + e, a)e-'^(^+^'^)'*da = / (x + e)(4-x-e 
Jo 
X + e)(4 — X — e) 



+ 00 



da 



i + (4-x-.) -^^-y/^ ' 



the unique solution is zo(x + e, x) = 1. Else, (15.120 admits a solution in [0, 1[. Writing that 
g{x + e, x) — 1 = edig{x,x) + o(e) (since g{x,x) = 1), and expanding the integrand of ()5.12l) 
with respect to e gives: 



1=1 + 



(^g^g^^^ 2.)^2(4 _ ^^2(-^ _ g-a^ 



x(4 — x) I xa + 4 — 2x 



e + o(e) > da 



=1 + 



1 dig{x, x) (4x5 - 48x^ + 192x3 - 256x2) + 48x3 - IBSx^ + 204x - 4x^ + 59 



4 x(x — 4)2 (x^ — 4x — 1 

The bracket vanishes for dig{x,x) given in (I5.10|) . From 

hD\zo{x,x)k{x,h)dh = / D\zQ{x,x)k{x,h)dh + I h\h\D^iZo{x,x)k{x,h)dh 



+ o{e). 



h'^ [—dig{x,x)],k{x,h)dh — / h"^ [dig{x,x)],k{x,h)dh, 



we obtain ()5.1ip . ■ 

The graph of x — > dig{x,x) is drawn in Figure [6] (b). We can verify that the fitness gradient 
vanishes at x* ~ 2.88. For x < x* the fitness gradient is positive, implying that the size tends 
to increase to x*. Similarly, for x > x*, the fitness gradient is negative and evolution reduces 
the size to x*. 
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Figure 6: (a): The Age- structured Canonical Equation, obtained from the Trait Substitution 
Sequence process in the limit of small mutation, (b): Graph of the function x i— > —dig{x,x) 
defined in i5.10\) . (c): The function ip defined in ^5.11\) . 



5.2 Comparison with a penalized logistic population without age-structure 



We compare the previous example with a logistic population, whithout age-structure, where 
the birth rate is b{x) = x(4 — x) (senescence is absent) and the death rate is d(x,M) = 1/4 + 
0.001(4 — x)M. The large population approximation corresponds to the logistic equation 

-^(x, t) = [ x(4 ~ ^) ~ ^ j "^-(2;, t) — 0.001(4 — x)m'^{x, t), m{x, 0) = mo(x), 

which admits the following unique solution in C^(M+,M+): 
(x(4 -x)-l) mo(x)e(^(^-^)~i/^)* 



m{x, t) 



(x(4 -x)-\)+ 0.001(4 - x)mo(x)e(^(''-^)-V4)i 



fh{x) :- 



[x(4- 



0.001(4 -x) 

(5.13) 

The infinitesimal generator of the TSS process takes here the following explicit form: 



L^{x) 



{x + h) — 4>{x))px 



[x(4-x)-l/4]. 
0.001 



1 



1 



4(4 - x)(4 - X - /i)y x + h 



h 



k{x, h)dh 



and the Canonical Equation becomes: 



dx 



px 
0.001 



X 4x(4 — x)2 



h^k{x, h)dh. 



A mutant y can invade the monomorphic resident population of trait x at equilibrium if 

^— ^ ^ > 4^ y - - y 4 - X - - y - X > 

y(4 -y) 4 

^y G ]min(x,/2(x)),max(x,/2(x))[ 



where /2(x) = 4 — (l/4)/(4 — x). The same conclusion as in Section 15.1.21 holds. The trait 
X* = 7/2 is an ESS to which the TSS and the Canonical Equation converge. 



Let us comment on these results. First, we point out that the ESS is larger in absence of senes- 
cence. In this case, the size (15.130 of the population at equilibrium for a given trait x is 
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larger since the birth rate does not decrease with age. The stronger competition then favors 
large sizes in the trade-off between growth and reproduction. Another reason for this difference 
is that senescence more or less reduces the reproduction period to the beginning of life. Without 
senescence, the individual has more flexibility in its reproduction strategy. It may choose a size 
that ensures him a longer life and that allows him to give birth at more spaced intervals. 

Let us also notice that the model with senescence is similar to a penalized version of the logistic 
model of this section, with a stronger natural death rate 

b{x) = x(4 - x), d{x, (Z, 1)) = ^ + 0.001(4 - x){Z, 1), 
in the sense that these models lead to the same equilibrium sizes ()5.5p and ESS x* (15. 9p . 

5.3 Age-logistic population interacting through the Kisdi interaction kernel 

We replace the logistic death rate by an age-logistic death rate where individuals interact through 
Kisdi's interaction kernel ()2.5p . In this example, computation becomes rapidly intricated. 

The constant Rq is the same as in ()5.2p . Equation (12. lip defining fh(x,a) becomes here 

9'rh , , /I Cva <r^ \ ^ , ^ ,^ ^ 

— {x,a)=[- + -^MAm{x,a). (5.14) 

As we have seen before, this type of equation is called age-logistic since the death rate in 
monomorphic populations is proportional to the size of the population. Plugging the solutions 
of (j5.14p into the boundary equation of (j2.12|) gives for a non trivial equilibrium: 

1= / x(4-x)e 4« 2" ^+-^ da = xU-x)e32c.M, — M _ $ / —]] 

Jo V Ci^M^ \ \^\jCuM^)) 

(5.15) 

where $ is the distribution function of the standard Gaussian law. Since the integral term defines 
a continuous and strictly decreasing function of M^, there exists for the traits satisfying (|5.2p 
a unique solution Mx to (j5.15p . There is however no explicit expression, and we use numerical 
computation to obtain the approximation presented in Figure El 




12 3 4 

Trait 



Figure 7: A numerical resolution of Equation !i5.15\) . giving for each trait x an approximation 
of the size Mx ■ 
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Once Mx has been computed, it is possible to obtain an approximation of the unique solution 
of ()5.14p . (I5.15P given as a function of by: 



m{x, a) 



Mr\ CuMr 



4 V CvM^ 



exp 



a CvMx 



4 2(1 + 1/) 



(5.16) 



To simulate approximations of the Age-structured TSS-jump process, we have to compute the 
fitness ZQ{y,x) of a mutant y in the monomorphic resident population of trait x. This can be 
obtained by solving numerically an equation in which the term is replaced by its numerical 
approximation M^: 



1 = 2/(4 - y) exp (-\a + {z- l)y{A - y){l - e^) - ^l^lh^a' j da. 



(5.17) 



(b) 





12 3 4 

Resident trait 



(c) 



(d) 



Figure 8: (a): Probability of extinction as a function of the resident and mutant traits, (b): 
Associated PIP. The domain of invasibility is delimited by the first bissector and by a curve that 
is not the graph of a function any more, (c): Simulations of an Age- structured TSS path, (d): 
Approximation of the Age- structured Canonical Equation obtained in the limit of the TSS when 
mutations are small. 

Conversely to the example corresponding to the rates (j2.3|) . (|2.4p . the TSS process does not 
seem to converge in long time here. The twisted curve in the PIP (Figure [8]) shows indeed that 
we are in a case of mutual invasibility (see [111 [26l 116]). 
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6 Example 2 



In this last example, we investigate a model where the death rate has the separable multiplicative 
form commented in Section [2.11 



(6.1) 



d{x,a,Z)= / a{l + e~'')U{x,y)Z{dy,da), 



U being the Kisdi kernel defined in (j2.5p . The birth rate is here b{x) = x(4 — x). Here, it will 
be possible to carry an investigation of the branching phenomenon that is observed. 

6.1 Monomorphic equilibrium in large population 

In this example, the partial differential equation defined in Proposition 12.41 is given by 

Qui QtII 

—-{x,a,t) + ——{x,a,t) = —am{x,a,t) / {1 + e~'^)U{x, x)m{x, a,t)da, (6.2) 
ot oa Jo 

(6.3) 



(6.4) 



m(x,0,t)= / x{4: — x)m{x,a,t)da. 
Jo 

Proposition 6.1. Equations Ii6.2\) - ^U7S\) admit a unique nontrivial stationary solution: 



■m(x, a) 



TTX^{A - xf{l + u) 



2Cu 



TTx'^{4 — x^a? 



where ^ is the distribution function of the standard gaussian law. 
Proof. p.l3p becomes here 

^ „ 2 ^ r+oo 

m{x, a) = m{x, 0)6" ^daE{x) ^ - ^y^Eix) ^ ^ ^ e-")?7(x, x)m{x, a)da. 

Jo 

(6.5) 

Plugging (j6.5p in (j6.3p . we obtain 

f + 00 



1 



E(^)a^ , ^ x(4-x) / 27r , 7rx2(4-x)^ ,^ ^, 
x(4-x)e 2 da 4^ 1 = — -\ — ^ E{x) = ^ —. (6.6) 







E{x) 



Finally, from (|6.5p 



m(x, 0) 



Eix) 



E{xf/^{l + i^) 



E(x)a'' 



/(,+°°(l + e-")C/(x,x)e —da Cv^fl^ 



i +e2sw 1 - $ 



(6.7) 



We deduce the announced result from (16. 5|) . (16. 6p and ([HI 
In Appendix [B1 we show that the nontrivial equilibrium ()6.4p is asymptotically stable. 
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(a) 



(b) 




Figure 9: (a): Extinction probability as function of the resident and mutant traits, (b): PIP 
showing the regions of possible invasion, obtained by solving Ii6.11\) . (c): Simulation of the Age- 
structured TSS process, (d): we show the successive equilibrium age distributions, and right, the 
sequence of corresponding traits. 

6.2 Invasibility 

Invasion of the resident population of trait x at equilibrium by a mutant of trait y is possible if 
and only if 

/ y(4 — y)e 2 da>l, where E{y,x) = / {I + e^")U{y,x)fh{x,a)da. (6.8) 
Jo Jo 

The balance equation (|2.14p is 

1= / y(4-y)e da, (6.9) 



Jo 

and hence, the probability of extinction ZQ^y, x) is strictly less than 1 if and only if: 

E{y,x)<E{y) (6.10) 

Since E{y,x) = -E(2;)^i±^^^^^r^, (|6.1U|) becomes: 

22 



This inequality can be solved numerically and the PIP is given in Figure [9] (b) . 

In the case where invasion is possible, the fitness of the mutant of trait y in the resident popu- 
lation of trait X is the smallest solution in [0, 1] of ()3.8p which in this example becomes 

1 = / y(4 - yy{^-mi~y)a-E(y,x)a^l2^^ 

Jo 

The solution of (|6.12p cannot be obtained explicitly, but can be computed numerically for every 
X and y in [0, 4] as obtained in Figure [9l 

Conversely to Example 1, the TSS process does not converge to a limit, even if it seems to stay in 
the neighborhood of the evolutionary singularity x* ~ 3.2. Moreover, looking at the underlying 
microscopic process (Figured]), we can observe splitting of the population into smaller groups 
once the neighborhood of this point has been reached. To understand this fact, let us study 
more carefully the directional derivatives of the solution of ()6.12p . 

(a) (b) 




X X 



Figure 10: (a): Graph of x ^ dig{x,x) defined in ^6.13\) . (b): Graph of x ^ d22g{x,x). (c), 
(d): graphs of x ^ d'lig{x,x) on (0,4) and on [3.2,3.4]. 
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Proposition 6.2. Recall that g has been defined in Lemma \3.5[ Then, for all x S (0,4), 

1.4(x + 1.4)(x-3.2) 

dig{x,x) = ^ 6.13 

x(4 — x) 

^2 , , 5.8(x2 + 1.8a; + I.l)(x2-6.6x + 10.9) 

dii9{x,x) = V2 

a;^(4 — xy 

2 . s_ 5.8(x2 + 1.82; + 1.5)(x2-6.6x + ll.l) 

2:^(4 — x)^ 

Proof. In order to obtain these derivatives, we replace as previously y hy y + Ey and x hy x + Sx 
in ()6.12p . expand the right hand side in series with Maple in Ex and Ey and identify the terms 
according to their order. ■ 

We remark that the function x dig{x, x) vanishes on the unique point x* = 3.2 of the interval 
[0,4], which confirms the simulations. d1ig{x,x) remains negative on [0,4] and the function 
y 1-^ g{y,x) — 1 is hence locally concave at the neighborhood of {x,x), for each x G (0,4). In 
particular, x* corresponds to a local minimum of the fitness function. Moreover, we can check 
that d22gix* , X*) ~ —2.9 < df^g{x* , x*) ~ —0.3 which entails that x* is a branching point in 
the terminology of [HI [26l E]. These results can also be read on the PIP (Figure [9] (b)) which 
shows that the population can still be invaded once it has reached the neighborhood of the 
evolutionary singularity x* . 



7 Conclusion 

In this work, we have considered a population structured by traits and by a continuous age, 
varying with time. The ecology is approximated by a partial differential equation in both trait 
and age. Based on a measure- valued process approach and generalizing [UIB], we have obtained 
new approximations extending the Trait Substitution Sequence and the Canonical Equation of 
Adaptive Dynamics to trait and age-structured populations. To our knowledge, these equations 
have not been proposed so far in the biological literature. Our approach emphasizes that the 
object of interest here is the whole age-distribution of the monomorphic population at equilib- 
rium, and not only the equilibrium trait. 



Adding an age structure to trait structured population models opens the way to new prob- 
lems dealing with life history feature. In the examples, we have seen that age-structure does 
not seem to change the qualitative behavior of the evolutionary approximation but refines the 
computation of the evolutionary stable state. As in the case without age-structure, the choice 
of the competition kernel has a strong influence on the patterns that can be observed on the 
microscopic simulations, and on the nature of the evolutionary singularities. 
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Appendix 

A Proof of Theorem 13.11 and of Proposition 13.31 

The skeleton of the proof of Theorem 13. II is similar to the one in Champagnat [1] for populations 
without age-structure, but with additional difficulties. In |4J, the distance between a population 
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and its equilibrium is obtained by comparisons of the trait values and of the sizes. Here, age- 
distributions have to be taken into account. We need to compare the measure with the 
stationary solution of (j2.8p . Thus, we consider the space of finite measures AiF{X) which we 
embed with the weak convergence topology, induced by the Dudley metric (see Rachev [33] p79): 

^H, u e Mf{^), D{fi,u) = sup 

f l-Lip 

ll/lloo < 1 

To prove Theorem 13.11 we establish that Vt G VF C X mesurable, Ve > 

^ hm^ P is monomorphic of trait y € T, D{Z^^^^^ ,Cy) < e) = {Xt S T) , (Al) 

where (Xt)tgK_^ has been defined in Theorem 13. 1[ We use the exponential deviations and esti- 
mates of times of exit of domains established in [37] one the one hand, and comparisons with 
linear age-structured birth and death processes on the other hand. This last point is detailed 
through the main steps of the proof of Proposition 13.31 in the sequel. For a complete proof, we 
refer to [36] (Chapter 6). 

Let £ > and let us introduce the following stopping times: 

i?;^=inf{t>0 I D[l{,^^ix)Z^idx,da),^,,idx,da)^ >£^ 

= inf |t > I l{y}{x)Zj'{dx, da) > 
= mi\t>0 \ I {x)Zl'{dx, da) = 



Recall that the time of first mutation r and the time of return to a monomorphic state 9 have 
been defined in Proposition 13. 3i By large deviations results, and by comparison results, we can 
prove that (see [36], Section 6.2.4) 

3p > 0, 3no G N*, Vn > no, Ko,q",y ( — < ^ ^ > 1 - 2e. (A2) 

On the period [0, r A i?" AS^ A Sq] , the mutant population is of mass smaller than e and evolves 
in a resident population that is closed to S^xo- The extinction probability ZQ{y, x) that appears in 
Theorem 13.11 comes from comparison of the mutant population process with the following linear 
birth and death age-structured (non-renormalized) processes: for sufficiently large n G N* (such 
that l-UnP>l- e), Vt G [0,r A i?^ A 5^], V/ G e,,(]R+,M+), 



^(Z^^/> < j^l{y}{x)f{a)Z^{dx,da) < -(^''',/>, 



where: 



1. {zl''^)t^K^ is the birth and death process with birth rate bi{a) = (1 — £)b{y, a) and death 
rate di{a) = d{y,a) + J^^ U{{y,a), {xo,a))m{xQ,a)da + Ue, 

2. {Z^''^)t^K^ is the birth and death process with birth rate 62(0) = b{y,a) and death rate 
^2(0) = d{y,a) + /jj^ U{{y,a), {xq, a))m{xo, a)da - Ue. 

Let us define for i G {1, 2} : 

5^„ = inf{t>0 I {Zi'',l)>en], 5^ = inf{t>0 | 1) = o} . 

The probabilities that these processes reach the level en before getting extinct are given by the 
following lemma, proved at the end of the section: 
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Lemma A.l. Let e > 0, and let i £ {1,2}. Let us consider the linear birth and death age- 
structured process {Zl''^)t^M.+ introduced above and starting from ZQ^{da) = 6o{da). Let 
be a positive real sequence such that lim„^_|_oo in/ log n = +00. 

1. // /+°° bi{a)e- fo M'^)'i'-da < I: 

lim P (si < t„ A Si^) = 1 (A3) 

n—y+oo 

2. If /+°° bi{a)e- ^0 M'^)'i^da > 1; 

lim P {Si, <tnA Sl^) = 4{y, xo) (A4) 
lim P <tr,<Sl) = \- xo), (A5) 

where ZQ{y,XQ) is the smallest solution in [0,1] of the equation analogous to i f 3. 5|) . where 
the birth and death rates are replaced by bi{a) and di{a). □ 

The probability that the mutant population gets extinct before having reached the mass e, 
before the occurrence of a new mutation, and before the resident population deviates from its 
equilibrium is then lower bounded by: 



^"o,.o">. (^0 <r^R:^ S-) > P^„,,.,, 5- < ^ A 5^ < r A 



>Ku<iS,y ( -^o < A Si, -^<rAR^]>zo{y, xq) - Ce, (A6) 

by ()A2|1 . by Lemma fA.ll and by showing that when e — > 0, Zo(y, x) -zo(y, x). The continuity of 
the extinction probability when the birth and death rates are perturbed by e is obtained thanks 
to the implicit function theorem (see Lemma 6.2.5 of ^S|)- Similarly, the probability that the 
mutant population reaches the mass e before extinction, before the occurrence of another mutant 
or before the deviation of the resident population to its equilibrium is lower bounded by: 



^xo,qS,y i^e < T A R^' A Sq) > P"o_gj_a;o ^5"" < A , 



< r A i?" 

nu„ 



>Ko,,E,^, (si <^aSI^<tArA>1- zo{y, xo) - Ce. (A7) 

If 2;o(y,xo) = 1, we obtain by (|A6p that: 

Pxo,,J,y (-^o < T A i?^ A S^) =PSo,,„",, < T A R^,, V, = xo, D{Z^, f^J < > 1 _ Ce. (A8) 

Since this is valid for every e > 0, Proposition 13.31 is proved in this case. If ZQ{y,XQ) g]0, 1[, 
we show by following the proof of Lemma 3 in [4J that once the mutant population has reached 
the mass e, it replaces the resident population with probability one. Indeed, the microscopic 
process in this case follows its deterministic dimorphic approximation. Since zo(?/,xo) < 1, we 
necessarily have (Proposition 13. 4( ): 

/ 6(y,a)exp(— / d(y, a, xo)da j da > 1, and (A9) 

by the assumption of non coexistence in the long term of two traits (Ass. 12. 6p : 

j 6(xo,a)exp ( — / d{xo,a,y)da \ da < 1. (AlO) 
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Then, the deterministic approximation ^ converges to 6y{dx)m{y,a)da. Its neighborhood is 
reached by the microscopic process in finite time with a probabihty that tends to 1 when n 
+00. When this happens, we can approximate the dynamics of the resident population by 
comparing it with a hnear birth and death age-structured process with birth rate b{xQ,a) and 
death rate d{xo,a,y) as we did for the mutant population after its introduction in the system. 
Because of (lAlOp . these linear branching processes can be chosen sub-critical and we can then 
show that the resident population gets extinct with a probability that tends to 1 when n — > +00. 
This gives us that: 

JimFl^^gj^^y ie<T, V = y)>l- zo{y,xo) - Ce. 



-00 



Since ()A6P gave us that: 

lim „n {e<T, V = xq) >zo{y, xq) - Ce, 
and since P^^.g^.j^ {0 < t, V = xq) + ^xo,q^,y ^ ^ = y) = 1> Proposition ESI is proved. 
Al Proof of Lemma I A. II 

For the linear birth and death processes, there is no accumulation of birth and death events and 



lim tn A Sin = +00, P - p.s. 

?i— >+oo 



By dominated convergence, the left hand side of (|A3p and (|A4p converges to P (5*0 < +00) , 
which solves an equation similar to ()3.5p . The result is then obtained from Proposition 13.41 
Let us now consider (IA5D: 



P {Sin <tn< S'q) =P {Sin < tn and S'q = +oo) + P (S'*^ < t„ < 5^ < +00) . (All) 

The second term of (lAlip is upper-bounded by P {tn < Sq < +00) , which tends to when 
n — > +00 by the choice of t„. Let us now turn to the first term. Under the assumptions of Point 
2, there exists a unique A > such that: 

Let Yi and W be the number of children and lifelength of an individual with birth rate 6j(a) and 
death rate di{a). From the sufficient conditions of Doney [H], when t —>■ +00, {e~^'^{Zl'^ , l))t£R+ 
converges almost surely and in mean square to a strictly positive random variable on {^q = +00} 

ifE{Y^) < +00 and E (^(^J^ e-^''bi{a)daj log (^J^ e~^%i{a)daj^ < +00. Recall Assumptions 
12.11 and let e > be sufficiently small so that d2 remains bounded below by d > 0. Since the 
density of W is di{w) exp{fQ di{a)da) and since conditionally to W, Yi is a Poisson random 
variable with parameter bi{a)da, we have 

/■+00 / /-w \ f+00 

E {Y^) = J i^J bi{a)daj di{w)e- -^o d.(o)da < + eU + UM^) J we-^'"dw < +00. 

The function x E M* 1-^ xlog(|x|) tends to when x — > 0. On M+, this function first decreases 
from to —e~^ and then increases. Thus: 



- < E 

e 



^ e-^%^{a)da^ log (^J^ e-^%i{a)da^ ^ < E (max(0, bW \og{bW))) 

r+oo 

<b{d + eU + UM^) / w \ln{bw)\ e-^'"dw < +00. 
Jo 
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Doney's sufficient conditions are then satisfied. On {Sq = +00} we thus have 

lim l^gi^ = A > 0. (A12) 

Let us consider n > 1/e, so that log(en) > 0. Since lim„^_|_oo Sl^ = +00 almost surely, we have 
by (|A12p that almost surely on {Sq = +00} 

hrn ^> hrn 'r =^> »• (^13) 

Then: 

lim P (Si^ < tn, si = +00) = hm P ( r^f^ < = +00 

=P(5^ = +oo) = l-4(y,x), (A14) 
as by choice of tn, lim„^_|_oo tn/ (Jog en) = +00. 

B Stability of the nontrivial equilibrium in Example 2 



Let us prove that the nontrivial equilibrium ()6.4p of Equations (I6.2p - ()6.3p is asymptotically 
stable. In view of Theorem 4.12 in Webb |^38j (p207), it is sufficient to study the eigenvalues 
of the operator B associated to the linearization of these equations in the neighborhood of this 
equilibrium. This leads us to study existence of non trivial solutions for the following equations: 

^(a) = - (aE{x) + X) u{a) - am{x, a)E{x), (Bl) 



da 

n(0) = / x(4 - x)u{a)da, (B2) 







-00 



E{x)= U{x,x){l + e-'^)u{a)da. (B3) 

Jo 

Equation ()Bip implies that: 

u{a) =e~^(^-)4-Aa J^^(o) - m(x, O)^(x) ae^'^do^ . (B4) 



Plugging this in Equation (|B2p : 

/+OO E( ) / f \ 

x(4 - x)e 1^-^" l^n(O) - m{x, 0)E{x) J ae^^da j da. (B5) 

If A = 0, then (1B5P and the balance condition (16. 9p lead us to: 

= m{x,0)E{x)x{4 — x) / — e 2 c?a, 

Jo 2 

which is never satisfied, since the right hand side is strictly positive. Hence A = is not an 
eigenvalue of B. 
If A / 0, JBS]) gives 

m(x, O)E(x) r+°° e"^^ (e"^" + Aa - l) da 
A2/o e-^(l-e-^«)da 
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If A > 0, then ()B6p has a sign opposite to E{x) and we deduce from (jB4p that so has u{a). 
This contradicts ()B3p unless both members of the equation are zero. This is excluded since u is 
an eigenfunction and hence is not constant equal to zero. Hence there is no eigenvalue on the 
nonnegative real axis. 

Let us now consider A = Ai + iX2 with Ai > 0. Since /q^°° ae~^^^^°'^ /'^da = 1/E{x), 
u(0) ^^(^'O)-^(^) 



and hence (IB4D becomes 



, . _E(x\^fh{x,Qi)E{x) \ 1 e-^'^ I 

^ ^ ^ i?(:r)/+°°e-T(l-e-^-)da^ 



Using ([BTD in ([BS]) 

0= / — — m(x,0)(l +e-'')e- 2 







1 e^'*''^ 
-a + --- \da-X=:A{X,x). (B8) 

We wish to check that for any given x in [0,4], the complex function A A(A,x) does not 




Figure Bl: The Jordan curves used to study the zeros of A i— > A(A, x) 

have roots in the half plane of complex numbers with nonnegative real part. Let x £ [0, 4] 
be fixed. When |A| +oo, the integral in (jBSp converges to a finite value. Hence, the roots of 
A(A, x) in "P"*", if they exist, belong necessarily to a compact set. In order to show that there are 
no roots on compact sets of "P"^, we use the argument principle (see Henrici [18], Section 4.10). 
The function A ^ A(A, x) is analytic on "P"*" \ {0}. Let F : r E [0, 1] i— > F(r) be a positively 
oriented Jordan curve of \ {0}. Assume that there are J zeros with multiplicities 

(mj)jgp J] contained in F. Then 

g = i / = h U„ = 

where „(A(r,.),0):=t5(^fe = M_jM, (bio) 

ZVr ZTT 
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where (/)(r) is a continuous version of the argument of A(r(r),x). Our purpose is to show that 
for chosen Jordan curves in \ {0}, the right member of (IBToI) is zero. This will entail that no 
zeros are contained in the chosen curves. We consider the Jordan curves drawn in Figure IB1[ 
Since the computations can not been carried explicitly, we compute numerically n(A(r,x),0) 
by following the algorithm proposed by Henrici ([18j, end of Section 4.6). We then let x vary 
along a grid between and 4. The numerical results tell us that the variation of the argument 
— i?!)(0))/(27r) remains zero for every x G [0,4]. We have thus checked numerically that 
there is no complex solution A of (|B8p with nonnegative real part, and the nontrivial equilibrium 
()6.4p is asymptotically stable. 
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